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ABSTRACT 

The spallation phenomenon was studied through numer- 
ical analysis using a coupled Lagrangian particle track- 
ing code and a hypersonic aerothermodynamics compu- 
tational fluid dynamics solver. The results show that car- 
bon emission from spalled particles results in a signif- 
icant modification of the gas composition of the post- 
shock layer. Preliminary results from a test-campaign at 
the NASA Langley HYMETS facility are presented. Us- 
ing an automated image processing of high-speed images, 
two-dimensional velocity vectors of the spalled particles 
were calculated. In a 30 second test at 100 W/cm 2 of 
cold- wall heat-flux, more than 1300 particles were de- 
tected, with an average velocity of 102 m/s, and most fre- 
quent observed velocity of 60 m/s. 

Key words: ablation; spallation; arc-jet; thermal protec- 
tion system. 

1. INTRODUCTION 

Low density ablators are one of the preferred class of 
material for protecting space vehicles during atmospheric 
entry. Through various physical processes [7, 14] such as 
near-surface oxidation, pyrolysis chemical reaction, va- 
porization, and other decomposition processes, these ma- 
terials reduce the amount of heat being conducted to the 
surface of the vehicle. These beneficial processes usu- 
ally fall under the umbrella term “ablation”. The spal- 
lation phenomenon [6, 11] also belongs in this category; 
however, unlike the other processes involved in ablation, 
spallation is generally regarded as undesirable. 
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Spallation is typically described as the ejection of solid 
particles off the surface of the material. Those parti- 
cle are then convected away in the flow field. The ex- 
act cause of this phenomenon is unclear, although it is 
likely to be caused by mechanical erosion of the fibrous 
surface, and weakened by oxidative decomposition, com- 
bined with the shear stresses introduced by the fluid flow. 
Soot formation from the pyrolysis gas is also a potential 
cause of spalled particles. 

The aerodynamic heating rates of re-entry vehicles may 
be affected by spallation through various mechanisms. 
For instance, the solid particles ejected from the surface 
can modify the near surface chemical processes and flow 
field. Mechanical erosion can lead to an irregular surface, 
and in turn, to an augmented surface roughness. Such a 
modification of the surface, as well as the presence of par- 
ticles in the boundary layer, could potentially trigger tran- 
sition from laminar to turbulent flow. Finally, mechan- 
ical erosion results in an accelerated surface recession, 
thereby reducing the effectiveness of the thermal protec- 
tion system (TPS). Material response codes, such as the 
ones used by NASA to design TPS, take into account this 
accelerated failure mechanism through introduction of an 
empirically calculated parameter [8]. Therefore, there is 
a need to better understand and model the spallation phe- 
nomenon, including its source and impact on TPS perfor- 
mance. 

To further investigate spallation effects, a Lagrangian par- 
ticle tracking code was developed. This code computes 
the dynamics of spalled particles as they travel through a 
flow field. The code uses a finite-rate chemistry model to 
predict the chemical interactions of the particles with the 
flow field. This particle tracking code can then be dynam- 
ically coupled to a computational fluid dynamics (CFD) 
solver which models the time-accurate hypersonic flow 
field around ablative samples. A preliminary analysis 
performed using this coupled approach is presented here. 

To support the spallation modeling effort, an experi- 
mental investigation was performed at the NASA Lan- 
gley Hypersonic Materials Environmental Test System 
(HYMETS) facility [4]. HYMETS is a 400 kW arcjet 



that can simulate hypersonic flight and Earth or Martian 
entry conditions. High speed imagery was used to capture 
spalled particles ejecting from FiberForm®, and PICA 
test samples at target heat flux conditions of 100, 200, and 
400 W/cm 2 . From this data, measurements of particle ve- 
locity and trajectory were obtained. An initial assessment 
of the data for one FiberForm® sample is presented here. 

2. NUMERICAL SIMULATIONS 
2.1. Numerical models 

The flow field solutions used to predict the trajectories 
of spalled particles following ejection from the surface 
were obtained using the hypersonic CFD code KATS 
- Kentucky Aerothermodynamics and Thermal-response 
Solver [16, 17]. KATS is a three-dimensional laminar 
Navier-Stokes solver that computes flow fields in thermo- 
chemical non-equilibrium in the continuum regime. The 
translational and rotational energy modes of the species 
of the gas are described by a single temperature T tr , while 
the vibrational and electronic energy modes, as well as 
electron translational energy mode, are characterized by 
another single temperature T ve . The individual species 
viscosity model used to compute the diffusive fluxes was 
obtained using the Blottner curve fits [1]. Euken’s rela- 
tions [13] were used to account for species thermal con- 
ductivity. The mixture transport properties are approx- 
imated using Wilke’s semi-empirical mixing rule. For 
reacting flows, a standard finite-rate chemical kinetics 
model was used. 

The spallation code used for this preliminary analysis 
models the dynamics of individual particles ejected from 
the surface of a test articles, using a steady-state flow- 
field solution as an initial condition. The code uses the 
flow field information to determine the trajectory of the 
particle. The following assumptions are used to derive 
the governing equations: 

1 . The spalled particles are spherical 

2. The spalled particles are made of graphite 

3. Heat and mass transfer rates are uniform over the 
surface of the particle 

4. Thermophysical properties of the particle are as- 
sumed to be uniform over the entire particle 

Graphite is chosen since charred PICA is considered to be 
carbon, and that thermochemical properties of graphite 
are available in the literature. A Lagrangian formulation 
is employed to compute the projected path of the spalled 
particle. The corresponding governing equations are 

-rhc/CO ~ ^C/CN - TTlC / sub\ 

f d 

Qconv H - Pdrag Qrad T Qrxn J 

(1) 

In the left hand side, m p , u p and E p are respectively the 
mass, the velocity vector and the energy of the spalled 
particle. In the right hand side, F D , q conv , p dr ag , Qrad 


Table 1. Surface chemistry model used to calculate the 
erosion of the spalled particles in the flow field 
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Table 2. Free stream and wall properties of the Mach 5 
argon flow 
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and q rxn are, respectively, the drag force, the convective 
heating, the power loss caused by drag, the re-radiative 
heating and the chemical energy term. 

As they travel in the flow field, the particles chemically 
react with the flow. As the particles lose mass, their radius 
is reduced. The gas/particle interaction is modeled using 
a finite-rate rate chemistry model that accounts for three 
types of non-homogeneous reactions: oxidation, nitrida- 
tion, and sublimation [3]. The three types are respectively 
represented in Eq. 1 by the mass source terms riic/eo > 
rile /cn and rile /sub- The complete reaction mechanism 
is provided in Table 1 . 

This spallation model and coupling strategy has been ex- 
tensively verified; more details are given in Ref. [3] 

2.2. Preliminary results 

The test case considered here consists of a Mach 5 high 
enthalpy flow over an axi-symmetric solid sample with a 
radius of 3.30 cm, and an axial length of 0.82 cm. These 
test conditions are representative of what is observed in 
the HYMETS facility [2, 4]. 

First, a parametric study was performed to evaluate the 
effects of the ejection conditions (ejection velocity, parti- 
cle size and ejection angle) on the trajectory of the parti- 
cle [3]. For the flow field, a Mach 5 argon flow was used, 
with flow conditions listed in Table 2. In the table, T w 
refers to the temperature at the surface of the sample. Ar- 
gon was chosen as it simplifies the flow field solution, and 
allow the particle to travel the flow field without reacting. 
As can be seen in Fig. 1(a) to Fig. 1(c), changing the ejec- 
tion conditions considerably changes the trajectory of the 
particle. Therefore, if trajectories are measured experi- 
mentally, they could provide valuable insights about the 
ejection conditions. 

Next, to allow evaluation of the effects of a spalled par- 
ticle interacting with a flow field upstream of a test arti- 
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(a) Particle radius variation (v p = 25 
m/s, 6 = 0°) 
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Figure 1. Particle trajectories as a result of the study for the ejection parameters of spalled particle in a Mach 5 argon 
flow 


cle, the spallation code was coupled with KATS. Using a 
time-accurate simulation, the carbon species emitted by 
the spalled particle traveling in the flow field was trans- 
ferred to the CFD code through mass source terms. For 
this specific case, a Mach 5 air flow field is used, with 
conditions listed in Table 3. In this table, Yi represents 
the mass faction of individual species i. For this simula- 
tion, a 30 micron carbon particle was ejected normal to 
the surface, 5 mm from the stagnation point, at a veloc- 
ity of 90 m/s. The simulation takes place over a physical 
time of 0.7 ms, and the results are presented in Fig. 2. 
In Fig. 2(a), the entire trajectory of the particle is seen 
over the temperature field. It can be seen in Fig. 2(b) to 
Fig. 2(d) that the particle ejects a significant amount of 
carbon species, and that this carbon diffuses around the 
particles in all directions before being convected away. It 
is clear from this preliminary analysis that these carbon 
species would react with the incoming flow, and that ad- 
ditional species such as CN could be created. These ob- 
servations reinforce the hypothesis that spalled particles 
could be responsible for a previously unexplained strong 
spectroscopic signature observed in some arc-jet test data 
upstream of the sample [5, 9, 10, 15]. 

3. EXPERIMENT 

3.1. Facility and instrumentation 

A series of preliminary test was performed at the 
HYMETS facility at the NASA Langley Research Cen- 
ter to assess the feasibility of identifying and tracking 
spalled particles ejected from TPS materials. As opposed 
to large-scale arc-jet wind tunnels, HYMETS can be op- 
erated by a minimal staff, can be run for a long period, 
and has essentially no downtime between runs. Hence, 


it is ideal for exploratory studies such as those of the 
present research effort. Moreover, the facility is equipped 
with numerous optical ports which allow the possibility 
of complex diagnostics and image capturing [4] . 

As previously noted, the purpose of the present test- 
campaign was to confirm the presence of spalled particles 
for different test articles throughout a range of flow con- 
ditions. During these tests, the facility was equipped with 
the following: 

• four high speed cameras for particle detection at dif- 
ferent view angles, two of them linked for stereo- 
scopic measurements 

• two spectrometers (VUV and N-IR) 

• an infrared camera 

• a two color single-point pyrometer for surface tem- 
perature measurements 

• three intrusive probes for flow calibration 

• three thermocouples at the back face of each sample 

A total of ten samples of PICA [12] and FiberForm® 
were tested in air plasma, under three cold- wall heat flux 
conditions (100, 200 and 400 W/cm 2 ). All samples had 
a diameter of 3.30 cm, and a height of 2.54 cm, which 
included a 0.38 cm Li-2200 collar for some samples. In 
addition, three graphite, and three Li-2200 samples were 
tested under the same heat fluxes. The following sec- 
tion present preliminary measurements on spalling par- 
ticle with high speed imaging, leaving a thorough doc- 
umentation of the test campaign measurements to future 
publications. 

Fig. 3 is a composite image compiled from all frames dur- 
ing a 30 s test of FiberForm® at 100 W/cm 2 . The images 


Table 3. Free stream temperature and wall properties of the Mach 5 air flow 
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(a) Trajectory and temperature 
field 


(b) Carbon density, t =0.17 
ms 


(c) Carbon density, t =0.35 
ms 


(d) Carbon density, t =0.53 
ms 


Figure 2. Temperature field , trajectory and density of the carbon species emitted from a 30 micron spalled particle ejected 
at 90 m/s in a Mach 5 airflow. 


used to construct this figure were obtained using a 12-bit 
monochrome pco.dimax HD high-speed CMOS camera, 
cropped down to 672x608 pixels from the 1920x1080 
full sensor array in order to reduce file size. As can 
be seen, numerous spalled particles are observed being 
ejected from the samples. The same observation was 
made at flows with lower cold- wall heat flux conditions. 
Spallation was systematically observed during all tests of 
FiberForm® and PICA samples. The graphite and Li- 
2200 samples however, showed no evidence of spallation. 


3.2. Particle tracking 

To provide a quantitative evaluation of the spalled parti- 
cle behavior, data was analyzed from an IDT MotionXtra 
N-4S1, 1016x1016, 10-bit mono camera located orthog- 
onally to the sample. The test case presented here is for a 
FiberForm® sample subjected to a 200 W/cm 2 heat flux 
for 30 s. The camera was set to double-exposure mode 
with the first exposure of 30 ps, followed 100 ns later by 
a second exposure of approximately 470 ps. Due to lim- 
itations in the on-board memory of the camera, both ex- 
posures were acquired at a sampling frequency of only 15 
Hz. This ensured that images were acquired for the entire 
duration of the 30 s placement of the sample within the 
arc-jet. As a result, 946 images were acquired over 31.5 
s, with 473 images acquired at each exposure duration. 



Figure 3. Integrated image of a FiberForm® samples 
subjected to a 400 W/cm 2 heat flux. The image is a re- 
construction formed by combining the brightest pixel of 
approximately 900 images taken during a 30-second ex- 
periment. These images were acquired at a frame rate of 
30 Hz with an exposure time of 25 /is 


It should be note that only the first exposure was ana- 
lyzed, with analysis of the second exposure left for future 
work. A single, unprocessed image is shown in Fig. 4(a). 
A small particle trace can be seen in the top of the im- 






(a) Original image (b) Enhanced image (c) Processed image 

Figure 4. FiberForm® sample subjected to a 200 W/cm 2 heat flux; after processing the image, six spalled particles are 
identified. 


age, and shadows of others at the bottom. However, by 
re-scaling the pixel intensities, as done in Fig. 4(b), multi- 
ple spalled particles are clearly observed, with five being 
identified in this single 30 ps exposure. 

To analyze the image, a processing scheme was devel- 
oped to isolate and enhance the particles. The first step in 
this process consists of subtracting the background infor- 
mation obtained by averaging the five consecutive images 
centered on the image of interest. This subtraction did not 
completely remove all of the radiated light emitted from 
the front of the sample, and therefore pixels within this re- 
gion were discarded. Next, using a pixel intensity thresh- 
old, the remaining pixels were identified and flagged as 
belonging to spalled particles. Finally, an algorithm was 
implemented which identified groups of neighboring pix- 
els within a 6x6 pixel moving window were identified 
and each separate group marked as a “spalled particle”. 
Groups containing ten or fewer pixels were identified as 
“false positives”. This process was fully automated using 
MATLAB® . The result is shown in Fig. 4(c). 

Using this process, more than 1300 particles were iden- 
tified over the 400 images obtained during a single run. 
By using the exposure length and the pixel grouping ge- 
ometry, two-dimensional velocity vectors in the plane of 
the camera can be estimated. A compilation of the veloc- 
ity vectors detected during the run is presented in Fig. 5. 
With this approach, the velocities are projected in the 
plane captured by the camera, and they are thus subjected 
to a two-dimensional approximation. As expected, the 
velocity around the centerline of the sample appear to be 
smaller since the velocity components in the z direction 
is not measured. 

From these results it is also possible to evaluate the nor- 
malized distribution of velocity for all the spalled par- 
ticles that were identified. Fig. 6 plots the distribution 
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Figure 5. Two-dimensional velocity vectors ( in the cam- 
era plane) of the detected spalled particles estimated by 
exposure length and pixel grouping geometry. Red line 
shows the location of the test article. 


of the particle with respect to their calculated xy- plane 
velocity. The probability is calculated by counting the 
number of particles which have a velocity within a cer- 
tain range, and then normalizing over the total number 
of particles. The average velocity component of particles 
in the two-dimensional image plane was estimated to be 
102 m/s, However the distribution of the particle veloc- 
ity was highly skewed, with the most frequent (probable) 
velocity closer to 60 m/s. 





Figure 6. Magnitude of the velocity distribution of the de- 
tected spalled particles. The magnitude is only calculated 
using two measured components, in the camera plane. 

4. CONCLUSION AND OUTLOOK 

This paper presents the recent advances in research re- 
lated to the ablative phenomenon - spallation. To inves- 
tigate the effects of spallation on the flow field, a nu- 
merical code was developed. The results showed that a 
spalled particle ejected in a reacting flow would release 
carbon species in the flow field. This carbon is in suf- 
ficient quantity that it is likely to form CN species, and 
therefore explain the CN emission detected by spectro- 
scopic experiments. 

High speed imaging measurements allowed observation 
of multiple spalled particles during an arc-jet test run on 
PICA and FiberForm® samples. An automated image 
analysis procedure allowed detection of as many as 1300 
single particles being ejected from the surface during a 
30 s test run. 

These results, although still preliminary, are very encour- 
aging. Particle tracking based on stereoscopic imaging 
will enable improvements in the determination of par- 
ticle velocities and trajectories in the three-dimensional 
space. Using the numerical method combined with the 
images obtained during the test campaign, it will be pos- 
sible to gather more information about the trajectory and 
initial state of the spalled particles. This will undoubtedly 
lead to a better understanding on how these particles are 
formed, and what process is responsible for their high- 
speed ejection. 
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